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Abstract 

In standard Regge theory, the contribution of the tripe-pomeron amphtude to 
the t = differential cross section for single diffraction dissociation has the form 
(fa/dM^dt\t=o ~ (s/so)^["^(°)"^V(A^^)"'^^°''> where apit) is the pomeron trajec- 
tory. For ajp(O) > 1, this form, which is based on factorization, does not scale with 
energy. From an analysis of pp and pp data from fixed target to collider energies, 
we find that such scaling actually holds, signaling a breakdown of factorization. 
Phenomenologically, this result can be obtained from a scaling law in diffraction, 
which is embedded in the hypothesis of pomeron flux renormalization introduced 
to unitarize the triple pomeron amplitude. 

t Fellow of CNPq, Brazil. 
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1 Introduction 



In Regge theory, the high energy behavior of hadronic cross sections is dominated by 
pomeron exchange fl], ^. For a simple pomeron pole, the pp elastic, total, and single 
diffractive (SD) cross sections can be written as 
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where apit) = a]p{0) +a't 



\+a't is the pomeron Regge trajectory, Pippp(t) is the 



coupling of the pomeron to the proton, g{t) is the triple-pomeron coupling, s' = is 
the IP —p center of mass energy squared, ^ = l—xp = s'/s = M"^ j s is the fraction of the 
momentum of the proton carried by the pomeron, and Sq is an energy scale parameter, 
which is assumed throughout this paper to be 1 GeV^ unless appearing explicitely. 

In analogy with Eq. ^, the term in brackets in (|^) is identified as the IP — p total 
cross section. 
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where we have used g{t) = g{0), since it was found experimentally that g{t) is indepen- 
dent of t [|1|. The remaining factor in (^, namely 



fip/p{^,t) = 



IGtt 



(5) 



where K = /5|,pp(0)/167r, is interpreted as the "pomeron flux". Thus, pp diffraction 
dissociation can be viewed as a process in which pomerons emitted by one of the protons 
interact with the other proton ||^. 

The function F{t) represents the proton form factor, which is obtained from elastic 
scattering. At small t, F'^[t) ~ e^'®* [Q. However, this simple exponential expression 
underestimates the cross section at large t. Donnachie and LandshofF proposed that 
the appropriate form factor for pp elastic and diffractive scattering is the isoscalar form 
factor measured in electron-nucleon scattering, namely 
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(6) 



where m is the mass of the proton. When using this form factor, the pomeron flux is 
referred to as the Donnachie-Landshoff (DL) fiuxQ. Note that at small-t Ff{t) can be ap- 
proximated with an exponential expression whose slope parameter, b{t) = ^ InFf (t), is 

^ The factor K in the DL flux is Kdl = (S/JiPqq )^/47r^, where /3]pqq is the pomeron-quark couphng. 
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4.6 GeV^^ at t ~ —0.04 GeV^, consistent with the slope obtained from elastic scattering 
at small t. 

As we discussed in a previous paper the ~ dependence of o"t(s) violates 
the unitarity based Froissart bound, which states that the total cross section cannot 
rise faster than ~ In^ s. Unitarity is also violated by the s-dependence of the ratio 
o'ei/o'T ~ s^, which eventually exceeds the black disc bound of one half (cXe; < ^ctt), as 
well as by the s-dependence of the integrated diffractive cross section, which increases 
with s as ~ s^^ and therefore grows faster than the total cross section. 

For both the elastic and total cross sections, unitarization can be achieved by eikon- 
alizing the elastic amplitude 0, 0], which takes into account rescattering effects. At- 
tempts to introduce rescattering in the diffractive amplitude by including cuts ||, Q or 
by eikonalization ^ have met with moderate success. Through such efforts, however, it 
has become clear that these "shadowing effects" or "screening corrections" affect mainly 
the normalization of the diffractive amplitude, leaving the form of the dependence 
almost unchanged. This feature is clearly present in the data, as demonstrated by the 



CDF Collaboration |T0[ in comparing their measured diffractive differential pp cross 
sections at a/s =546 and 1800 GeV with pp cross sections at ^/s = 20 GeV. 

Motivated by these theoretical results and by the trend observed in the data, a 
phenomenological approach to unitarization of the diffractive amplitude was proposed 
based on "renormalizing" the pomeron flux by requiring its integral over all ^ and t to 
saturate at unity. Such a renormalization, which corresponds to a maximum of one 
pomeron per proton, leads to interpreting the pomeron flux as a probability density 
simply describing the ^ and t distributions of the exchanged pomeron in a diffractive 
process (see details in section 3). 

In this paper, we show that the hypothesis of flux renormalization provides a good 
description not only of the s-dependence of the total integrated SD cross section, as 
was already shown in Q], but also of the differential (or ^) and t distributions. 
Specifically, we show that for > 5 GeV^ (above the resonance region) and ^ < 0.1 
(the coherence region [|l|), all available data for p + p{j>) — > X + p{j>) at small t can be 
described by a renormalized triple-pomeron exchange amplitude, plus a non-diffractive 
contribution from a "reggeized" pion exchange amplitude, whose normalization is fixed 
at the value determined from charge exchange experiments, pp — *■ Xn. A good fit to 
the data is obtained using only one free parameter, namely the triple-pomeron coupling, 
^7(0). 

We also show that the t = cross section at small ^ displays a striking scaling be- 
havior, namely cPa/dM'^dt\t=o ~ C/(M^)"-^^°'', where the coefficient C is s-independent 
over six orders of magnitude. In contrast, the ~ s^^ dependence expected from the stan- 
dard triple-pomeron amplitude represents an increase of a factor of 6.5 between y/s = 20 
and 1800 GeV. This scaling behavior is predicted by the renormalized flux hypothesis 
and provides a stringent and successful test of its validity. 

In section 2 we present and discuss the data we use in this paper; in section 3 
we describe our phenomenological approach in fitting the data using the pomeron flux 
renormalization and pion exchange models; in section 4 we present the results of our fits 
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to data; in section 5 we present the case for a scaling law in diffraction; and in section 
6 we make some concluding remarks on factorization and scaling in soft diffraction. 



2 Data 



The data we use are from fixed target pp experiments 0|, from ISR pp experi- 



ments [|T3|, from SppS Collider pp experiments and from pp experiments at the 
Tevatron Collider [0, |l6l. Below, we discuss some aspects of the Tevatron Collider data 



reported by CDF ^ . 



2.1 The CDF data 

The CDF Collaboration reported ||T0| differential cross sections for pp — >■ pX in the 
region of ^ < 0.15 and \t\ <~ 0.15 GeV^ at -y/i = 546 and 1800 GeV. The experiment 
was performed by measuring the momentum of the recoil antiproton using a roman pot 
spectrometer. No tables containing data points are given in the CDF publication. The 
data are presented in two figures (figures 13 and 14 in flOl), which are reproduced here 
as Figs, la and lb, where the number of events is shown as a function of xp [xp = 1—0 
of the recoil antiproton, rather than as a funcion of the antiproton momentum (as was 
done in [|1^). The histogram superimposed on the data in each figure is the CDF fit to 
the data generated by a Monte Carlo (MC) simulation. As an input to the simulation, 
the following formula was used: 
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The first term in this equation is the triple-pomeron term of Eq. ^ The second term 
was introduced to account for the non-diffractive background. A connection to Regge 
theory may be made by observing that 7 = 1 (0) corresponds to pion (reggeon) exchange 
with a Regge trajectory of intercept a(0) = (0.5) (see section 3). The factor of | does 
not appear in reference [10| and is introduced here to account for the fact that we refer 
to the cross section for pp — > pX and do not include that for pp —>■ pX, as was done 
by CDF. The CDF MC simulation took into account the detector acceptance and the 
momentum resolution of the spectrometer. The slope of the pomeron trajectory, a', was 
kept fixed at the value a' = 0.25 GeV~^. The values of the remaining parameters, as 
determined from the CDF fits to the data, are listed in Table 1, where we include the 
values for the momentum resolution, ctq, at ^/s = 546 and 1800 GeV. 



2.1.1 Acceptance corrected .^-distributions 

Using the information provided in the CDF publication, we mapped Figs, la and lb into 
Figs. 2a and 2b, respectively, in which the data are corrected for detector acceptance. 
The acceptance was obtained from Fig. 2 of reference [|l^. The results are presented as 
cross sections, rather than as events, versus x. The normalization was determined by 
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comparing the data points with the CDF MC fits. The number of events corresponding 
to each x-bin of the MC histograms in Figs, la and lb was converted to an absolute cross 
section by convoluting the analytic CDF formula for the differential cross section with the 
t-acceptance function and with the Gaussian ^-resolution function using a normalization 
that reproduces the MC histogram. The curves in the new figures represent Eq. |^ 
convoluted with a Gaussian resolution function of ^, whose width was determined from 
the momentum resolution of the spectrometer at each energy. Specifically, these curves 



are calculated using the expression [|T7| 
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where is given by Eq. |^ (with ^ and g{^', ^) is the Gaussian resolution function 
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As seen in Figs. 2a and 2b, expression (|]) provides an excellent fit to the acceptance- 
corrected differential cross sections, including the unphysical region of negative values. 
Thus, once the detector experimental resolution is accounted for, the low ^ (or equiva- 
lently, the low M^) cross section is completely compatible with that expected from ex- 
trapolating the cross section from the region of 0.95 < xp < 0.99 (0.05 > ^ > 0.01) into 
the resolution dominated very low-^ region using the triple-pomeron differential cross 
section shape. This behavior rules out the hypothesis of low ^ (low M^) suppression 



suggested by some authors [|19| |2T|, |2^ . 



2.1.2 Cross sections at t 



-0.05 GeV^ 



As can be seen in Fig. 2 of reference []10[, the CDF data in the triple-pomeron dominated 
region of ^ < 0.05 are concentrated at low t-values, namely |t| <~ 0.1 (0.2) GeV^ 
for = 546 (1800) GeV. Therefore, direct comparison of the CDF data with other 
experiments should be made for t-values within these regions of t. Since the CDF paper 
does not report ^-distributions at a fixed value of t in the form of a table, we extracted 
such a table for t = —0.05 GeV^ from the information given in the CDF paper. The 
value of t = —0.05 GeV^ was chosen in order to allow direct comparison of the CDF 
data with the data of reference |rT|, for which ^-distributions have been published in 
table form for t = -0.05 GeV^ and = 14 and 20 GeV (see Tables 2 and 3). The 
t = —0.05 GeV^ CDF points were evaluated from the data in Figs. 2a and 2b, which 
represent cross sections integrated over t, by scaling the cross section at each point in C, 
by the ratio 

d'^a/d^dt\t=-om 
= d^^ ' 

which was calculated using Eq. ^. Figures 3a and 3b display the t = —0.05 GeV^ data 
points grouped into ^-bins of approximately equal width in a logarithmic scale. Fig- 
ures 3c and 3d display in a linear ^-scale the data for ^ < 0.01, including the unphysical 
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region of negative ^-values. The horizontal "error bars" represent bin widths. The val- 
ues of the points plotted in Figs. 3a-d are listed in Tables 4 and 5. The solid (dashed) 
curves in the figures represent the CDF fits without (with) the convoluted ^-resolution 
function, calculated using Eq. (Eq. |). For = 546 (1800) GeV, the effect of the 
detector resolution becomes important for ^ < 0.005 (0.003). Immediately below these 
values, the data lie higher than the extrapolation of the solid-line fits from the larger 
^-values (see Figs. 3a,b). This effect is completely accounted for by the smearing effect 
of the ^-resolution, which also accounts for the values of the cross sections in the un- 
physical negative ^-regions, as seen in Figs. 3c and 3d (for exact numerical comparisons 
see Tables 4, 5). The effect of the resolution on the measured cross sections is quite 
substantial at low ^ and therefore must be taken into consideration when comparing the 
low-^ CDF data with predictions of unitarization models based on low-^ suppression of 
the diffractive cross section [jl9|, EOl pTl . 



2.1.3 t-dependence 

We now return to the question of the slopes bo of the t-distributions (see Table 1). The- 
oretically, the value of bo for pp pX should be the same at all energies and equal to 
one half of the corresponding value for pp —* pp (see Eqs. and Experimentally, 
lbo,ei = 4.6 GeV-2 [|. The best-fit CDF slope values are bo = 7.7 ± 0.6 (4.2 ± 0.5) 
GeV-2 for ,/s = 546 (1800) GeV. The 1800 GeV value is close to 4.6, within error, but 
the 546 GeV slope is significantly larger than 4.6 GeV~^. The discrepancy between the 
slope value measured by CDF at y/s = 546 GeV and the expected value of bo = 4.6 
GeV~^ may be explained by the very short t-range of the experimental measurement. In 
the region of low ^, where pomeron exchange is dominant, the detector had reasonable 
acceptance only within the region 0.03 < |t| < 0.1 GeV^. Thus, the slope could not be 
measured accurately. The quoted error in the measured slope is the standard deviation 
calculated keeping all other parameters fixed at their best-fit values. The large corre- 
lation coefficients between the error of the CDF best-fit parameter bo and other fit 
parameters indicate that a good fit to the data within the t-region of the measurement 
could have been obtained with a different value of bo, and correspondingly different value 
of the other parameters, subject to the constraint that the integrated cross section over 
the t-range of the measurement remain the same. Since t = —0.05 GeV^ corresponds 
approximately to the cross-section-weighted mean value of t in the region 0.03 < t < 0.1, 
the value of the differential cross sections at t = —0.05 GeV^ is insensitive to a change 
in bo. 



2.1.4 Total diffractive cross sections 

At y/s =546 (1800) GeV, the total integrated cross section within the region > t > — oo 
and (1.5 GeV^) /s < ^ < 0.05 calculated using Eq. |^ (multiplied by a factor of 2 to include 
the cross section for pp Xp) is 7.28 (8.73) mb. 
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3 Phenomenological approach 



In the framework of Regge theory |23|, the cross section for pp — > pX in the region of 
large s/M\ can be expressed as a sum of contributions from exchanges of reggeons i, j 
and k (see Fig. 4), 

- '^G.,,{t)[^] [^] cos[Ht)-m] (11) 



dMldt s^f-^j^ '^^^'{Mj^J \ So 



with 



Gijkit) — f^ippi't) Pjppi't) Pkppi^) gijkit), (12) 

where ai{t) = ai{0) + a[t is a reggeon trajectory, [3ipp is the reggeon couphng to the 
proton, Qijk is the "triple- reggeon" coupling and (piit) is a phase factor determined by 
the signature factor, rfiit) = -)_ e-^'^^K*)^ where C = ±1 is the signature of the exchange. 
The signature factors have been expressed as r]i{t) = rf^it) e**^'*^*) with the moduli rf^it) 
absorbed into the (3 (t) parameters in (|12|) . For pp —>■ pX reggeons i and j must have the 
same signature, so that — (pj{t) = |[aj(t) — aj(t)]. As mentioned in section 1, the 
energy scale Sq is not determined by the theory and is usually set to 1 GeV^. The lack 
of theoretical input about the value of Sq introduces an uncertainty in the pomeron flux 
normalization, which is resolved in the renormalized pomeron flux model (see discussion 
below) . 

Table 6 displays the s and ^, or M|-, dependence of the contributions to the pp pX 
cross section at t = coming from various combinations of exchanged reggeons. Three 
Regge trajectories are considered: the pomeron, IP, with ajp(O) = 1 + e, the reggeon, 
R, with aR = 0.5, and the pion, tt, with = 0. In fitting elastic and total cross 
sections, Covolan, Montanha and Goulianos use two reggeon trajectories, one for the 
f/a family with a//a(0) = 0.68 and the other for the p/u family with ap/^{0) = 0.46; 
Donnachie and Landshoff use one "effective" trajectory with a'j(^{0) = 0.55. For 
simplicity and clarity in presentation, we consider in Table 6 one reggeon trajectory 
with a/j(0) = 0.5. The terms IPIPIP (triple-pomeron) and IPIPR correspond to the 
picture in which pomerons emitted by one proton interact with the other proton 
to produce the diffractive event. The last row in Table 6 shows the predictions of the 
renormalized pomeron flux model, which is discussed below. 



3.1 Standard approach 

The standard approach to diffraction is to perform a simultaneous fit to the pp —* pX 
differential cross sections of all available data at all energies using Eq. |Tl|, which is 
based on factorization. In such a fit, the only free parameters are the tripple-reggeon 
couplings, gijkit). The reggeon trajectories and the couplings Pit) are determined from 
the elastic and total pp cross sections [^, and the coupling P^^ppit) is obtained from the 
coupling P-irnpit) measured in the charge exchange reaction pp nX; using isotopic spin 
symmetry, p^ppit) = ^P^mpit)- 
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3.1.1 Factorization 



Equation [Tl| is based on factorization. A "global" fit of this form to all available data was 



performed by R. D. Field and G. C. Fox in 1974 [^. However, the data available at that 
time could not constrain the fit well enough to test the triple-reggeon phenomenology, 
let alone determine the triple-reggeon couplings. By 1983, with more data available from 
Fermilab fixed target and ISR experiments [|rT|-|n|, good fits to the small-t differential 
pp — > pX cross sections were obtained using the empirical expression |jl| 

The first term in ( 1T3|) has the ^-dependence of the IPIPIP amplitude with ajp(O) = 1 
(e = 0) and the second term has the ^-dependence of the nnlP amplitude. Note that a 
reggeon-exchange contribution, RRP, with q;_r(0) = 0.5, would have a fiat ^-dependence. 
At the relatively low values of -/s of the Fermilab fixed target and ISR experiments, 
the ^-range was not large enough for the 6-slope to be sensitive to the variation with ^ 
expected from Eq. 0, namely b = Bq — 2a' In^, or to distinguish between a 1/^ and a 
l/^^"*"^ dependence in the first term of (|13|) and thereby establish the now well known 
deviation of ajp(O) from unity. Nevertheless, the prominent ~ 1/.^ behavior of the cross 
section at low showed IPIPIP dominance and left little room for contributions from 
other terms, as for example from a IPIPR term with its sharper ~ i/^i+o.s dependence 
on ^. This is illustrated by the fits of Eq. ^ to the very precise data for pd — > Xd shown 
in Fig. 5. The data [^, ^ are from the experiment of the USA-USSR Collaboration 
at Fermilab using an internal gas-jet target operated with deuterium. The values of the 
cross sections at t = —0.035 GeV^ plotted in Fig. 5 were obtained either directly from 
the published Tables or by extrapolation from their published values at t = —0.05 
GeV^ using the measured slope of the t-distribution. The two sets of data were 
normalized to the average value of the cross section within the ^-region common to both 
sets of data. Figures 5a and 5b show fits using a 1/^ and a l/S}^" dependence (with 
e = 0.104 0), respectively. Both fits are in good agreement with the data. 

In summary, the agreement of the Fermilab fixed target and ISR experimental results 
with the empirical expression (p!3D, which is inspired by the factorization based standard 
tripple-reggeon phenomenology, shows that: 

• At small ^ the cross section is dominated by the IPIPIP amplitude (~ 1/0- 

• At larger ^ there is an additional contribution, which has the form of the -ktiIP 
amplitude (~ 0- 



3.1.2 Breakdown of factorization 

In 1994, when CDF published the diffractive pp pX cross sections at ^/s = 546 
and 1800 GeV |T3, the supercritical pomeron trajectory with ajp(O) > 1 was already 
well established by fits to total hadronic cross sections Therefore, CDF made fits 
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using Eq. 0, which includes two terms: the IPIPIP amphtude (first term) and a non- 
diffractive contribution parameterized as (Pand/dC,dt = IS,"'e^^. The form of the latter 
was inspired by the empirical expression (|T3p, and the parameter 7 was introduced to 
gj^^eciwe/?/ incorporate possible contributions both from TrvriP (7 = 1) and RRP (7 = 0) 
amplitudes, as discussed in section 2.1. 

Three important results from the CDF fits to the data should be emphasized: 

• Only the IPIPIP term and a non-diffractive contribution are required by the fits. 
An upper limit of 15% was set on a possible contribution of a IPIPR term to the 
total diffractive cross section at a/s = 546 GeV. From this result, we derive the 
following limit for the ratio, i?, of the coefficients (jjjfc(O) of the IPIPR/ IPIPIP 
terms: 

'-^ipipipxy) 



This limit is 5% of the value of R used in the fit by Erham and Schlein [19, 20| 
(see also comments in 



The parameter e was determined for the first time from the ^-distribution of single 
diffraction dissociation and was compared to the e obtained from the s-dependence 
of the total pp cross section |23|. The CDF results are: 



e (from total cross section) = 0.112 ± 0.013 (15) 
e (from rfd/rf^; ys = 546 GeV) =0.121 ±0.011 (16) 
e (from dd/d^; = 1800 GeV) =0.103 ±0.017 (17) 

The values obtained from the da/d^ distributions are, within the quoted uncertain- 
ties, consistent with the value determined from the rise of the total cross section, 
as would be expected for pomeron pole dominance. The weighted average of all 
three values is [Q 

e = 0.115 ±0.008. (18) 

Using the relation ^ = M^/s, the IPIPIP (first) term in Eq. ^ can be written in 
terms of as (see also Eq. |lT]) 



d^a G{0) . , .s s 



dM^dt 2 ^ ' (M2) 



s/s,Y . j^^^ e(^°+2"' (19) 



where in standard Regge theory 6 = 2e. Treating 5 as a free parameter and 
performing a simultaneous fit to the diffractive cross sections, agd, at i/i = 20 [pT 
546 and 1800 GeV, CDF obtained S = 0.030 ± 0.016. 



The last result indicates a breakdown of factorization. The observed slower than ~ (s/sq)^^ 
increase of the diffractive cross section with energy is necessary to preserve unitarity 
and was predicted in 1986 0| by calculations including shadowing effects from multi- 
ple pomeron exchanges. More recent work based on eikonalization of the diffractive 
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amplitude or on the inclusion of cuts shows that shadowing can produce substan- 
tial damping of the s-dependence of the cross section but has no appreciable effect on 
the M^-dependence. These predictions are in general agreement with the conclusions 
reached by the CDF fits to data. However, the damping predicted by the eikonaliza- 
tion model is not sufficient to account for the observed s-dependence of the total single 
diffraction cross section (see Fig. 6); the predictions of the model based on cuts are in 
better agreement with the data [^. 



3.2 Renormalized pomeron flux approach 
3.2.1 Triple-pomeron renormalization 

The CDF measurements showed that, just like at Fermilab fixed target and ISR energies, 
the shape of the small (small ^) behavior of the diffractive cross section at the 



Tevatron Collider is described well by the IPIPIP amplitude displayed in Eq. |I9|. The 
total diffractive cross section, obtained by integrating Eq. |T9| over all t and over 
from M^-^ = 1.5 GeV^ to M^^^ = 0.1s, increases with s as ~ {s/sqY. For 5 = 2e, 
which is the value for simple pole exchange, agd would increase faster than the total 
pp cross section, which varies as ~ {s/sqY, leading to violation of unitarity. With the 
experimentally determined value of 5 ~ 0, the diffractive cross section remains safely 
below the total cross section as s increases, preserving unitarity. 

As we have seen in the previous section, introducing shadowing corrections can 
dampen the increase of the diffractive cross section with s and thereby achieve the 
desired unitarization while preserving the M^-dependence of the IPIPIP amplitude, as 
required by the data. However, the shadowing models do not account completely for the 
s-dependence of the data, and the two models mentioned above do not predict the same 
amount of s-damping of the cross section. In addition, these models are very cumber- 
some to use in calculations of single diffraction, double diffraction and double-pomeron 
exchange processes. 

The calculational difficulties of unitarity corrections in the standard approach are 
overcome in the "pomeron flux renormalization" approach proposed by Goulianos 0. 
The renormalized flux approach is based on a hypothesis, rather than on an actual 
calculation of unitarity corrections, and therefore can be stated as an axiom: 



The pomeron flux integrated over all phase space saturates at unity. 



The standard pomeron flux is displayed in Eq. ^ Using F'^{t) = e^"*, the integral of 
the standard flux, 

N{s)= / / fjp/,{i,t)didt, (20) 



is given by 



N{s) = K—[Ei{r-2ehiUin)-E,{r-2t\YiUa.)i (21) 
2a 
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where Ei{x) is the exponential integral function^, r = boe/a', ^min = Mq/s = 1.5/ s is 
the effective diffractive threshold, and ^max = 0.1 0. 

The renormalized pomeron ffux, f^i^, t), can now be expressed in terms of the stan- 
dard flux, fip/p{C,,t), as follows^: 

t) - I f^/^^^^'^ ifiv(.)<i 

The renormalized IPIPIP contribution to the differential cross section is given by 
or, in terms of M^, by 



dW^dt N{s) (M2)i+2^ 



■ {My. (24) 



In the energy interval of ^/s = 20 to 2000 GeV, the standard flux integral varies as 
~ s^*^ (see Fig. 7). Thus, flux renormalization approximately cancels the s-dependence 
in Eq. ^ resulting in a slowly rising total diffractive cross section. Asymptotically, as 
s CX3, the renormalized total diffractive cross section reaches a constant value: 

hm af,(s) = hm ^ = 2a!^^e^/^ (25) 

s^oo '^"^ ' s^oo N(s) 

The s-dependence of the integral of expression (^) over all t and ^ < 0.05, multi- 
plied by a factor of 2 to account for both pp — * pX and pp —>■ Xp, is compared with 
experimental data for asd{^ < 0.05) in Fig. 6 (from [Q). In view of the systematic 
uncertainties in the normalization of different sets of data, which are of (9(10%), the 
agreement is excellent. 

3.2.2 Pion exchange contribution 



The form of the empirical expression ([131) suggests that at high ^ the dominant non- 
IPIPIP concontribution to the cross section comes from pion exchange. In Regge theory, 
the pion exchange contribution has the form 

§i = U(M-.-'is(), (26) 

where fn/p{^,t) is the pion flux and cr'^^(s^) the irp total cross section. 



E^{x) = 7 + Ina; + where 7 =0.57721. . . (Euler's constant) 

For a detailed dis( 
which N{s) = 1 see [Q 



n=l 

For a detailed discussion of the role of the scale parameter sq in determining the value of s for 
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In the "reggeized" one-pion-exchange model |^^, the pion flux is given by 



Jn/p[^,t) = — 



\t\ 



Ait Att {t ~ mlY 



(27) 



where g^^p/ATi ^ 14.6 p4| is the on mass-shell coupling, a.„(t) = 0.9 1 is the pion trajec- 



tory, and G\{t) is a form factor introduced to account for off mass-shell corrections. For 
Giit) we use the expression (see [p9| and references therein) 



G{t) 



2.3 - ml 
2.3 -t 



{2i 



Since the exchanged pions are not far off-mass shell, we use the on-shell vrp total 
cross section [Q, 



7r ' p 



TT p\ 



10.83 (sO°-^°^ + 27.13 (sO 



-0.32 



(29) 



3.2.3 A one parameter fit to diffraction 

Motivated by the success of the empirical expression (|13]) in describing the Fermilab fixed 
target and ISR data, and by the similarity between this expression and the CDF fits 
to data at Tevatron energies, we have performed a simultaneous fit to single diffraction 
differential cross sections at all energies using the formula 

|| = /^(e, t) ■ a^^^isO + /./p(e, t) ■ a-^isO , (30) 

in which the first term is the renormalized triple-pomeron amplitude, Eq. and the 
second term is the pion exchange contribution, Eq. ^ Results from our fit, in which 
only the triple-pomeron coupling, gpipip, is treated as a free parameter, are presented 
in the next section. 



4 Results 

In this section, we present the results of fits performed to experimental data using Eq. |30| , 
which has two contributions: a renormalized triple-pomeron amplitude and a reggeized 
pion exchange term. 



4.1 Differential cross sections 



The experimental ^-distributions are usually distorted in the low-^ region by the reso- 
lution in the measurement of the momentum of the recoil p{p). We therefore check first 
how well Eq. |3y reproduces the shapes of the differential cross sections of the pp data of 
E396 JlTl at ^/s=U and 20 GeV and of the pp data of CDF ^ at yi=546 and 1800 
GeV in the regions of ^ not affected by detector resolution. Figure 8 shows the cross 
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sections dPa^d/didt at t = -0.05 GeV^ for E396 and CDF (data from Tables 2, 3, 4 



and 5). The solid lines represent the best fit to the data at each energy using Eq. |3G 
with the normalizations of the triple-pomeron and pion exchange contributions treated 
as free parameters. The quality of these fits indicates that no reggeon terms other than 
the triple-pomeron and pion exchange terms are needed to describe the shapes of the 
differential ^-distributions. 

Figures 9-10 show the result of a simultaneous fit (solid lines) to the t = —0.05 GeV^ 
data of E396 and CDF using Eq. |3y with only the triple-pomeron coupling as a free 
parameter. The overall normalization of the data was allowed to vary within ±10% to 
account for possible systematic effects in the experimental measurements. The shift in 
the normalization of the data at each energy that resulted in the best fit is given in each 
plot. In Fig. 10 the individual contributions of the triple-pomeron and pion exchange 
terms are shown by dashed curves. The fit had a = 1-0 per degree of freedom. 

The parameters used in the fit are aip{t) = 1.104 + 0.25t and /3jppp(0) = 6.57 GeV^^ 
(4.1 mba) for the triple-pomeron term, and those given in section 3.2.2 for the pion ex- 
change term. The fit yielded a triple-pomeron coupling gjppp = 1.0 GeV~^ (0.62 mb^), 
which corresponds to = 2.6 mb; using the Fi(t) form factor yields gipipip = 1.1 GeV" 
(0.69 mb^) and = 2.8 mb. 



Figure 11 shows a fit of Eq. ^to ISR data [|T3] of da/d^ versus ^ at fixed t. In this 
fit, the experimental ^-resolution was taken into account by convoluting Eq. ^ with 
the Gaussian resolution function using ctq = 0.003. The parameters used in Eq. [3^ 
were those derived above. The overall normalization of the data has an experimental 
systematic uncertainty of 15% 

4.2 Total difFractive cross sections 

In Fig. 12, we compare experimental results for the total diffractive cross section within 
< -t < -oo and ^ = M|./s < 0.05 with the cross section calculated from the triple- 
pomeron term of Eq. ^ (solid line) using the triple-pomeron couplig evaluated from our 
fit to the differential cross sections. Within this region of ^, the expected contribution 
of the pion exchange term is less than 2% at any given energy. The data points are from 
references [|TU], 0, 0, 0, |T|, |T|, |T§. 



There are two points that must be kept in mind in comparing data with theory: 

• Normalization of data sets 

The overall normalization uncertainty in each experiment is of (9(10%). 

• Corrections applied to data 

Deriving the total cross section from experimental data invariably involves ex- 
trapolations in t and ^ from the regions of the measurement to regions where no 
data exist. In making such extrapolations, certain assumptions are made about 
the shape of the t-distribution and/or the shape of the ^ distribution. Different 
experiments make different assumptions. For example, with the exception of the 



ISR experiments 13, 14], all measurements of the experiments listed here are at 
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very low-t. In these experiments, an exponential form factor of the form e^o^ is 
assumed for extrapolating into the high-t region. The (higher-t) ISR data show 
a clear deviation from exponential behavior and support the -Ff (t) form factor. 
Using -Ff (t) instead of e^o^ results in a larger total cross section by ~ 5 — 10%, 
depending on the value of s (smaller correction at higher s). The magnitude of the 
correction depends on the ^-region (and through ^ on s), since the t-distribution 
depends not only on the form factor but also on ^ throught the term e*^"^"'^"^^)*. 
Another source of error comes from the fact that the slope of the t-distribution 
is usually not measured accurately in experiments sensitive only to low-t. The 
discussion in section 2.1.3 of the CDF measurement at ^/s = 546 GeV illustrates 
this point. 

Table 7 presents the total diffractive cross sections corrected for the effects mentioned 



above. The ISR [T^, |T^ and SppS cross sections were left unchanged, since they were 



calculated taking into account the high-t behavior of the differential cross section. The 



cross sections of Refs. |T^, |T6| were multiplied by the ratio of the cross section calculated 
from the IP IP IP expression using the Fi{t) form factor to that calculated using the 
simple exponential form factor. Finally, the cross sections of Refs. [|T0|, |11| , for which the 
data are within a limited t-region and have no reliable slope parameters, were calculated 
as follows: in each case, we evaluated the integrated cross section within the t — ^ 
region of the experiment using the parameters determined by the experiment, and then 
recalculated this cross section using the formula of Eq. adjusting the normalization 
parameter D to obtain the same value for the integrated cross section over the same t — C, 
region; this formula was then integrated over the region < |t| < oo and 1.5/s < < 
0.05. The corrections to values derived directly from the published data are of (9(10%). 
In view of the systematic uncertainties in the normalization of the various data sets, as 
evidenced by the discrepancies among data from different experiments in overlapping 
s-regions. Fig. 12 shows excellent agreement between the experimental cross sections 
and the predictions of the one-parameter fit of Eq. ^ (using the Fi{t) form factor and 
do^P = 2.8 mb). 



5 A scaling law in diffraction 

The renormalization of the pomeron flux to its integral over all available phase space may 
be viewed as a scaling law in diffraction, which serves to unitarize the triple-pomeron 
amplitude at the expense of factorization. 

The breakdown of factorization is illustrated in Fig. 13, where cross sections are 
plotted as a function of ^ at fixed t for -^i = 14 and 20 GeV (\/Ii=17 GeV) and 
\/s2=1800 GeV. As noted by CDF [|l^, while the shapes of the da/d^ distributions as ^ 
decreases tend to the shape expected from triple pomeron dominance at both energies, 
the normalization of the S2 points is approximately a factor of {S2/S1Y = 2.6 lower than 
that of the si points, instead of being a factor of {S2/S1Y higher, as one would expect 
from factorization (see factor s' = in Eq. |^). 



14 



This particular way in which factorization breaks down imphes that the dPa / dMxdt\t=o 
distribution is approximately independent of s, and therefore scales with s, in contrast 
to the s^^ behavior expected from factorization. Figure 14 shows the differential cross 
sections as a function of M|. at t = -0.05 GeV^ yi=14, 20, 546 and 1800 GeV within 
^ regions not including the resonance region of Mx < 5 GeV^ (a/s=14 and 20 GeV) 
and not affected by the detector resolution (for i/s=546 and 1800 GeV). These cross 
sections are also shown in Fig. 15 for regions of ^ low enough not to be affected by the 
non-pomeron contribution. In Fig. 15, the data are compared with a straight line fit of 
the form da/dMx ~ l/Mj^^^", (solid line) and with the predictions based on factorization 
(dashed lines). Clearly, factorization breaks down in a way that gives rise to a scaling 
behavior. 

The scaling of the M'j^ distribution is a consequence of the pomeron fiux renormaliza- 
tion hypothesis, as pointed out in section 3.2.1. Figure 7 shows that the renormalization 
factor based on fiux scaling has an approximate s^^ dependence, which cancels the s^*^ 
dependence in da/dMx expected from factorization. An exact comparison between data 
and theory is made in Fig. 16, where data and predictions of Eq. |30| are shown for t=0. 
The t = data were obtained from the t = —0.05 GeV^ data shown in Fig. 14 by 
subtracting the pion exchange contribution at t = —0.05 GeV^ and calculating the t = 
cross section assuming a t-distribution given by F^^(t)e''~^"'^'^^^*. The excellent agreement 
between data and theory over six orders of magnitude justifies our viewing the pomeron 
fiux renormalization hypothesis as a scaling law in diffraction. 

6 Conclusions 

We have shown that experimental data on diffractive differential cross sections d^cr / d^dt 
for pp — * Xp and pp — > Xp at energies from ^/s =14 to 1800 GeV, as well as total 
diffractive cross sections (integrated over ^ and t), are described well by a renormalized 
triple-pomeron amplitude and a reggeized pion exchange contribution, whose normal- 
ization is kept fixed at the value determined from pp —>■ Xn. 

The renormalization of the triple-pomeron amplitude consists in dividing the pomeron 
fiux of the standard Regge-theory amplitude by its integral over all available phase space 
in ^ and t. Such a division provides an unambiguous normalization of the pomeron fiux, 
since the energy scale factor, sq, which is implicit in the definition of the pomeron pro- 
ton coupling /3jppp(0) that determines the normalization for the standard fiux, drops 
out. Thus, the renormalized pomeron fiux depends only on the value of ^min and on 
the pomeron trajectory, which is obtained from fits to elastic and total cross sections. 
Therefore the only free parameter in the renormalized triple-pomeron contribution to 
soft diffraction is the triple-pomeron coupling constant, gjpjpjp. From our fit to the data 
we obtained the value gpipip = 1.1 GeV~^. 

The scaling of the pomeron fiux to its integral represents a scaling law in diffraction, 
which unitarizes the diffractive amplitude at the expense of factorization. A spectacu- 
lar graphical representation of this scaling is provided by the experimental differential 
da/dMx\t=o distribution as a function of Mx for energies from y^=14 to 1800 GeV. 
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This distribution shows a clear ~ 1/(M|^)^+^ behavior, which is independent of s over 
six orders of magnitude, in agreement with expectations from the flux renormalization 
hypothesis and contrary to the ~ s^^ behavior expected from the standard theory based 
on factorization. 
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Table 1: CDF fit-parameters from reference [10]. 





x/i = 546 GeV 


y/s = 1800 GeV 


D 

bo 
a' 
e 
I 

7 
b' 

CTO 


3.53 ±0.35 

7.7 ±0.6 
0.25 ± 0.02 
0.121 ±0.011 

c-q7+498 
"J-J' -280 

0.71 ±0.22 
10.2 ± 1.5 
1.4 10-^ 


2.54 ±0.43 

4.2 ±0.5 
0.25 ± 0.02 

0.103 ±0.017 
162l^f 
0.1 ±0.16 

7.3 ±1.0 
8.9 10-^ 



Table 2: Differential cross sections at 14 Table 3: Differential cross sections at 20 
GeV and \t\ = 0.05 GeV^ [11]. GeV and \t\ = 0.05 GeV^ [11]. 





(Pa/d^dt 
(mb/GeV^) 






d^a/d^dt 
{mb/GeV^) 


0.0160 


282.0 ± 11.8 




0.0160 


233.2 ± 10.9 


0.0267 


145.6 ±9.3 




0.0267 


146.7 ±7.9 


0.0373 


112.0 ±8.4 




0.0373 


105.9 ±7.0 


0.0480 


100.0 ±7.7 




0.0480 


78.8 ±6.2 


0.0586 


85.8 ±7.3 




0.0586 


80.7 ±6.5 


0.0693 


79.7 ±7.1 




0.0693 


70.2 ±6.1 


0.0800 


69.1 ±7.6 




0.0800 


57.0 ±5.8 


0.0906 


65.4 ± 7.4 




0.0906 


62.5 ±6.6 


0.1013 


51.0 ±7.5 




0.1013 


68.6 ±7.0 
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Table 4: Differential cross sections at 546 GeV and \t\ — 0.05 GeV^ extracted from the 
CDF measurements [10] (see text for details). 





(Pa/d^dt 
(mb/GeV^) 


CDF Fit 

(mb/GeV^) 


Fit(8)Gauss 
{mb/GeV^) 


-0.0046 


21.18 ±4.73 




16.95 


-0.0027 


539.1 ± 40.3 




569.1 


-0.0009 


3534.4 ± 124.0 




3591.8 


0.0009 


4561.3 ± 135.1 


2568.1 


4561.3 


0.0027 


1682.9 ±73.8 


772.5 


1618.2 


0.0046 


563.0 ± 38.8 


443.4 


536.2 


0.0064 


300.9 ± 30.5 


308.7 


329.0 


0.0082 


226.5 ±25.7 


236.3 


244.8 


0.0100 


178.2 ±22.6 


191.5 


195.9 


0.0119 


136.5 ±20.6 


161.3 


163.9 


0.0146 


107.3 ± 13.9 


131.0 


132.3 


0.0192 


82.4 ± 11.3 


101.0 


101.6 


0.0256 


77.9 ± 11.7 


78.5 


78.8 


0.0348 


67.0 ± 11.6 


62.5 


62.6 


0.0458 


49.6 ±9.8 


53.3 


53.3 


0.0577 


54.4 ± 10.6 


48.6 


48.6 


0.0714 


41.2 ±7.8 


46.4 


46.4 


0.0870 


47.7 ±7.9 


45.9 


45.9 


0.109 


44.5 ± 7.1 


47.0 


47.0 
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Table 5: Differential cross sections at 1800 GeV and \t\ — 0.05 GeV^ extracted from the 
CDF measurements [10] (see text for details). 





(Pa/d^dt 
(mb/GeV^) 


CDF Fit 

(mb/GeV^) 


Fit(8)Gauss 
{mb/GeV^) 


-0.0022 


375.4 ± 48.0 




307.8 


-0.0011 


3419.4 ± 182.8 




3419.4 


0.0000 


8368.6 ± 278.9 




8368.6 


0.0011 


5646.9 ± 210.4 


1603.4 


5019.4 


0.0022 


1311.9 ±88.4 


776.6 


1311.9 


0.0033 


568.7 ±66.5 


513.7 


573.6 


0.0044 


403.4 ± 57.7 


386.1 


404.5 


0.0055 


319.6 ±52.9 


311.3 


320.0 


0.0072 


222.7 ±35.8 


243.7 


247.4 


0.0100 


196.7 ±33.8 


182.9 


184.2 


0.0139 


153.6 ±29.3 


140.1 


140.6 


0.0189 


106.7 ±22.1 


112.1 


112.3 


0.0250 


84.6 ± 18.8 


93.8 


93.9 


0.0322 


90.2 ± 18.7 


81.6 


81.7 


0.0422 


73.9 ± 13.7 


72.2 


72.3 


0.0555 


55.0 ±9.4 


65.3 


65.3 


0.0717 


69.9 ± 10.0 


60.8 


60.8 


0.0918 


57.6 ±7.2 


57.8 


57.8 


0.116 


55.4 ±6.5 


55.8 


55.8 
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Table 6: Triple-reggeon amplitudes for pp —* pX . 
ap(0) = 1 + e ajj(O) = 0.6 a,(0) = 

Amplitude (fa/d£,dt\t=o d^a/dM'^dt\t=o a*fo{s) 
{]P]P)P ~p5^ 
{]P]P)R ~pg4^ 

{RR}P ~ s'C ~ ~ 

{RR)R ~p|^ 

(iPi?)iP 

Renormalized: 

^(M^ -constant 
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Table 7: Total single diffraction cross sections for ^ < 0.05. The cross sections of the 
references marked with f were derived from the experimental data using the procedure 
outlined in the text. 



[GeV] 



c^sd [mb] 
(both sides) 



Ref. 



14 
20 
16.2 
17.6 
19.1 
23.8 
27.2 
23.4 
26.9 
30.5 
32.3 
35.2 
38.3 
23.3 
27.4 
32.4 
35.5 
38.5 
44.7 
53.7 
62.3 

546 
1800 

546 
1800 



3.94 ±0.20 
4.46 ±0.25 
4.87 ±0.08 
4.96 ±0.08 
4.94 ±0.08 
5.19 ±0.08 
5.42 ±0.09 

6.07 ±0.17 
6.05 ±0.22 
6.37 ±0.15 
6.32 ±0.22 
7.01 ±0.28 

6.08 ±0.29 
6.5 ±0.2 
6.3 ±0.2 
6.5 ±0.2 
7.5 ±0.5 
7.3 ±0.4 

7.3 ±0.3 
7.0 ±0.3 
7.5 ±0.3 

9.4 ±0.7 
8.46 ± 1.77 
8.34 ±0.36 
9.12 ±0.46 



Q 
H 

m 
m 
m 
m 
m 
m 
m 
m 
m 
m 
m 
m 
[0 







m 

m 
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CDF Single D iff r active Data 

(uncorrected for acceptance) 
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Figures la (top) and lb (bottom) 

CDF single diffraction data (uncorrected for acceptance): number of events versus xp] 
the solid line histograms are from a Monte Carlo simutation using formula (7). 
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CDF Single Diffractive Cross Sections 

(corrected for acceptance) 
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Figures 2a (top) and 2b (bottom) 

CDF cross sections da/dxp (integrated over t); the solid curves represent formula (7) 

and the dashed curves formula (8). 
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CDF Single Diffractive Cross Sections 

{da/di at t = -0.05 GeV^) 
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Figures 3a (top) and 3b (bottom) 

CDF cross sections da/d^dt at t = —0.05 GeV ^; the solid curves represent formula 

and the dashed curves formula (8). 
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CDF Single Diffractive Cross Sections 

{da/di at t = -0.05 GeV^) 
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Figures 3c (top) and 3d (bottom) 

CDF cross sections da/d^dt t — —0.05 GeV ^; the solid curves represent formula 

and the dashed curves formula (8). 
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Figure 4: Illustration of the triple-reggeon phenomenology. 



27 



Cross Sections for p -\- d ^ X -\- d 
• da/d^ at t = -0.035 GeV^) 
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Figures 5 (a,b): Cross sections for p + d — > X + d 
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Figures 6: The total single diffraction cross section for p{p) + p ^ p{p) + X vs i/i 
compared with the predictions of the renormalized pomeron flux model of Goulianos ^ 
(solid line) and of the model of Gotsman, Levin and Maor (dashed line, labeled 
GLM); the latter, which includes "screening corrections", is normalized to the average 
value of the two CDF measurements at ^/s = 546 and 1800 GeV. 
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Figure 7 : The integral of the standard pomeron flux for pp pX, N{s) of Eq. [2^ 
using -F^(t) = e^'^^, as a function of y/s (solid curve) is compared with a dependence 
~ s^^ (dashed curve). The horizontal solid line at N{s) = 1 represents the saturated 
renormalized flux. The flux integral calculated using in Eq. ^ the Fi (t) form factor of 
Eq. P can be approximated by the expression 0.41 s^*^. 
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Cross sections (Pagd/d^dt sd, t = 
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Figure 8 : Cross sections dPagd/d^dt for p + p{p) p{p) + X at t = —0.05 GeV^ and 
^/s = 14, 20, 546 and 1800 GeV. The solid lines represent the best fit to the data at 
each energy using two terms, the IPPIP and TTTriP amplitudes, with their normalizations 
treated as free parameters. 
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Figure 9 : Cross sections ^ ■ dPagd/d^dt for p + p{p) + X at t = —0.05 GeV^ 

and ^/s = 14, 20, 546 and 1800 GeV are compared with the results (sohd hues) of 
a simultaneous one parameter fit with a renormalized IPIPIP amplitude and a pion 
exchange contribution. To account for systematic uncertainties, the normalization of 
each data set was allowed to vary within ±10% of its nominal value; the parameter 
"N(data)" represents the shift in the data normalization for which the best fit was 
obtained. 
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Figure 10 : Cross sections (Pagd/d^dt for p + p{p) — > p{p) + X aX. t = —0.05 GeV^ 
and -y/s = 14, 20, 546 and 1800 GeV are compared with the results (sohd hues) of 
a simultaneous one parameter fit with a renormalized IPIPIP amplitude and a pion 
exchange contribution. The dashed lines represent the individual pomeron and pion 
contributions. To account for systematic uncertainties, the normalization of each data 
set was allowed to vary within ±10% of its nominal value; the parameter "N(data)" 
represents the shift in the data normalization for which the best fit was obtained. 
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Figure 11 : Cross sections (Pagd/ d^dt for pp pX measured at the ISR at various 
values of and t. The sohd hues are fits obtained using the renormahzed IPIPIP amph- 
tude and the pion exchange contribution, convoluted with the experimental ^-resolution, 
which dominates the shape of the distributions at small ^. The overall normalization of 
the data has a systematic uncertainty of 15% [13|. 
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Figures 12: Total single diffraction cross sections for p{p) + p — > p{p) + X versus 
^/s compared with triple-pomeron predictions based (a) on pomeron pole dominance in 
standard Regge theory (dashed line) and (b) on the renormalized pomeron flux model |Q 
(solid line). The cross sections were corrected for effects due to extrapolations in t, as 
discussed in the text. The errors shown are statistical; typical systematic uncertainties 
for each experiment are of (9(10%). 
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Figure 13 : Cross sections dPag^/d^dt for p + — >■ p{p) -\- X aA, t — —0.05 GeV^ and 
^/s = 14, 20 and 1800 GeV. The solid lines are the global one-parameter fit to the data 
presented in Fig. 10, and the dashed lines represent the renormalized triple-pomeron 
contribution. The dashed-dotted line is the standard flux triple-pomeron contribution 
at = 1800 GeV predicted from the data at = 14 and 20 GeV. 
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Figure 14 : Cross sections (Pagd/dM'^dt for p + p(p) +X ait — —0.05 GeV^ and 

y/s = 14, 20, 546 and 1800 GeV. 
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Figure 15 : Cross sections cPagd/ dM'^dt for p + p{p) p{p) + X at t = —0.05 GeV^ and 
^/s = 14, 20, 546 and 1800 GeV. At y^=14 and 20 GeV, the standard and renormalized 
flux fits have the same normahzation; at the higher energies, the renormahzed and 
standard flux predictions are shown by the solid and dashed hues, respectively. 
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Figure 16 : Cross sections d^agd/dM^dt for p+p{p) p{p)+X at t = and ^/s = 14, 20, 
546 and 1800 GeV, multiplied by N{s)/\/3'^ s'^''{P]pppg]p]p]p)/16n], are compared with 
the renormahzed flux prediction of 1/(M^)^+^. The dashed curves show the standard 
flux predictions. The t = data were obtained by extrapolation from their t — —0.05 
GeV^ values after subtracting the pion exchange contribution. 
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